Colloids in active fluids: Anomalous micro-rheology and negative drag 
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We simulate an experiment in which a colloidal probe is pulled through an active nematic fluid. 
We find that the drag on the particle is non-Stokesian (not proportional to its radius). Strikingly, 
a large enough particle in contractile fluid (such as an actomyosin gel) can show negative viscous 
drag in steady state: the particle moves in the opposite direction to the externally applied force. 
We explain this, and the qualitative trends seen in our simulations, in terms of the disruption of 
orientational order around the probe particle and the resulting modifications to the active stress. 
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Active particles take energy from their surroundings 
and convert this into mechanical work. Active fluids are 
suspensions of such particles in a Newtonian solvent, and 
they represent an interesting class of nonequilibrium soft 
matter [U [2j. To lowest order, an active fluid may be 
modelled as a collection of force dipoles, exerted by the 
active particles, which quite often have orientational or- 
der (e.g., as a nematic phase). These force dipoles are 
either contractile, when the forces are exerted "inwards" 
towards the centre of mass of each particle, or extensile, 
in the opposite case. Suspensions of bacteria such as 
E.coli are extensile fluids, while a dispersion of Chlamy- 
domonas (algae) is contractile, as is the actomyosin gel 
which constitutes the cytoskeleton of eukaryotic cells [l] . 

The continuous dipolar forcing present in active flu- 
ids profoundly affects their macroscopic properties. For 
instance, active nematics flow spontaneously in steady 
state in the absence of any external force, provided the 
activity level (dipole density) is high enough [3l|4]. Simu- 
lations [3l[5] show this flow to be chaotic, resembling the 
"bacterial turbulence" observed in concentrated films of 
B. suhtilis [6]. Furthermore, the bulk rheology of active 
fluids is strongly non-Newtonian |2l |7] . Theory predicts 
an increase (reduction) in the effective viscosity for con- 
tractile (extensile) fluids [2] . These predictions were con- 
firmed by simulations [H [9] , and also by experiments on 
Chlamydomonas [TO] and E. coli [TT]. 

The local flow properties of active fluids also differ from 
their passive counterparts. These can be addressed by 
monitoring the dynamics of a probe particle in a "micro- 
rheology" experiment. Such studies [12 find marked vi- 
olations of the fluctuation-dissipation theorem [TsHIl], 
which in near-equilibrium systems links the decay of ran- 
dom fluctuations to the linear force response. Local flow 
of active fluids is of strong biophysical relevance: for in- 
stance, the cytoskeleton in moving cells is mainly sub- 
jected to localized cues rather than global forces, and its 
response to these may be crucial to cell motility [16]. 

Here we address a very basic issue of the local flow re- 
sponse in active fluids. We ask: what happens when a 
passive particle of radius R is dragged through an active 
nematic with a force F7 This represents perhaps the sim- 



plest micro-rheological experiment possible on an active 
fluid. We find by simulations that this simple experiment 
should lead to some very interesting results. First, the 
drag coefficient ^ = F/v (with v a steady particle speed) 
exhibits a strongly nonlinear dependence on radius, in 
violation of Stokes' law. This is especially noticeable in 
the contractile case. We explain these violations in terms 
of the deformation of the orientational order in the ac- 
tive fluid around the probe particle, and present a simple 
scaling theory for the balance between active and viscous 
forces. Second, and strikingly, our theory and simulations 
show that in a contractile fluid the colloidal probe, if large 
enough, should move steadily in the direction opposite to 
the applied force. That is, we find a stable steady state 
with a negative drag coefficient. We also simulate a tran- 
sient (force reversal) experiment that probes further the 
remarkable physics in this regime. 

Numerical model: In the continuum limit, the hydro- 
dynamics of an active nematic fluid can be described 
by a set of continuum equations [TJ [2] that govern the 
time evolution of the velocity field and of a (trace- 
less, symmetric) tensor order parameter Qafd- The lat- 
ter describes the orientational order of the active parti- 
cles (whether bacteria, algae, or cytoskeletal filaments) 
which usually have a rod-like shape and are thus capable 
of nematic alignment [1 . Without activity, nematics are 
described by a Landau - de Gennes free-energy density 
^ = F{Qa^) + K{d^Qa^)y2, with 

(1) 

where indices denote Cartesian coordinates, summation 
over repeated indices is implied, 7 controls the magnitude 
of nematic order, and K is an elastic constant. 

The hydrodynamic equation for the evolution of the 
order parameter is: DtQa(3 = ^Ha(3^ with Dt a ma- 
terial derivative describing advection by the fluid ve- 
locity and rotation/stretch by flow gradients (see 
[3 ). The molecular field is Ha^s = —ST/SQa/s + 
/ ^Qa/s) and r is an inverse rotational fric- 
tion. The fluid velocity obeys daU^ = 0, and also 
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the Navier-Stokes equation, in which a passive Qa^- 
dependent thermodynamic stress enters [3^. The active 
force dipoles then create a further stress, 
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where ( is the activity parameter that sets the dipolar 
force density [2 . Note that ( < for contractile flu- 
ids and ( > for extensile ones. Within our hybrid 
numerical scheme, we solve the Navier-Stokes equation 
via lattice Boltzmann, and the equation for the order 
parameter via finite difference [17 . Periodic boundary 
conditions are deployed throughout. 

We introduce a spherical colloidal probe by a standard 
method of bounce-back on links [18 which provides a 
no-slip boundary condition for the velocity field at the 
solid surface. Order parameter variations create addi- 
tional elastic and active forces on the particle which are 
computed by integrating the total stress tensor over its 
surface. There we impose planar anchoring of the ne- 
matic director, which as usual is a headless unit vector 
ria oriented along the major principal axis of Qai3' Planar 
anchoring is imposed via an additional quadratic term in 
the free energy (see ^91 (20] for details). 

Below we give our results in simulation units [21]. To 
convert them into physical ones, relevant for instance to 
a contractile actomyosin solution, we can assume a value 
of K of 1.25 nN, and a rotational viscosity of 1.0 kPa/s. 
(These values hold for typical cytoskeletal gels [22].) In 
this way, the simulation units for force, length and time 
can be mapped onto 25 nN, 0.5 jam and 10 ms respec- 
tively. Note that the same equations also apply to an 
extensile bacterial fiuid, but the mapping to physical 
units in that case leads to very different values [3 . Our 
model also neglects motility which is important in bac- 
terial fiuids, where it naturally leads to density inhomo- 
geneities [3l [23] . 

Results: We now discuss the result of our micro- 
rheological simulation in which a colloid of radius R is 
pulled through an active fiuid, either contractile or ex- 
tensile. The external force F was directed either along 
or perpendicular to the far-field nematic director. For 
the values of R chosen, plots of the steady state veloc- 
ity V (measured along F) versus F (see Fig. la,b) show 
a well-defined linear regime at small external force. In 
all our simulations ( is kept small enough such that no 
spontaneous fiow arises in the absence of the probe par- 
ticle [25]. Figs. la,b show that the linear drag coefficient 
^0 = F/v\f^o increases with activity for a contractile 
fluid, and decreases for an extensile one (within the pa- 
rameter range that we explored) . Thus a contractile fluid 
opposes motion more strongly than its passive counter- 
part; an extensile fluid less. So far, this is in line with the 
respective increase (contractile) and decrease (extensile) 
of bulk fluid viscosities mentioned above ^4] . 




FIG. 1: (Color online) (a,b) Plots of v{F) curves for col- 
loids of fixed radius R — 11.3, at various ^ (see legend) in a 
contractile (a) and an extensile (b) active nematic, dragged 
respectively parallel and perpendicular to the far-field direc- 
tor. (c,d) Dependence of ?7eff, measured according to Eq. [s] 
for both (c) contractile and (d) extensile cases. Filled sym- 
bols refer to an external force parallel to the far-field director, 
while empty symbols are used for orthogonal drag. Different 
symbols (colours) refer to different activities (see legend in 
(a) and (b)). In (a) and (b), insets, we sketch the director 
field and active force directions for a colloid pulled (a) in a 
contractile active fluid along the far-field director, and (b) in 
an extensile fluid perpendicular to it. 



By analogy with the formula for Stokes drag on a probe 
of radius in a passive fluid of viscosity r^, we can deflne 
an effective viscosity for the probe motion as 



^eff(i^) =^0/(6^i?). 



(3) 



This is plotted as a function of R in Figs. lc,d after cor- 
rection for periodic boundary effects [26 . If Stokes' law 
did hold, the effective viscosity should be independent of 
i?, and we see that this is indeed the case for a passive ne- 
matic (C = 0, open and fllled circles in Fig. lc,d). When 
activity is switched on, this picture changes dramatically 
and the drag coefficient becomes strongly non-Stokesian. 
For a contractile ffuid, this anomaly is most apparent 
when pulling along the director fleld, where the effec- 
tive viscosity increases sharply with the radius. In the 
extensile case the most markedly anomalous response is 
instead obtained when the particle is pulled perpendicu- 
larly to the far-fleld director, and the effective viscosity 
this time drops with size. 

These results show that for active fluids microrheol- 
ogy experiments do not simply probe the fluid, but mea- 
sure a scale-dependent property of the probe and sur- 
rounding fluid in combination. Although microrheology 
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is thereby disconnected from bulk rheology, such exper- 
iments may offer essential insights into the physical re- 
sponse of (for instance) an actomyosin gel on the length 
scale relevant to transport of organelles or other subcel- 
lular objects [16 . 

We can qualitatively explain the anomalous drag 
^o(^, C) ill Fig. Ic via the following argument (sketched 
in Fig. la, inset). As the colloidal probe is pulled through 
the contractile fluid, it deforms the director field. With 
the pulling force oriented parallel to the far-field director, 
planar anchoring requires an elastic splay of the direc- 
tor locally. Because of the fluid flow created by particle 
displacement, this splay will be stronger in front of the 
moving particle than behind it. Through its effect on the 
active stress, the splay will result in a large force opposing 
the motion at the front, and in a smaller force favoring it 
at the rear. The net effect is to slow the particle down: 
^0 is increased. The same argument holds for pulling per- 
pendicular to the director, but here is much smaller. We 
attribute the different magnitude of the effect in the two 
directions to the fact that contractile activity enhances 
the splay response (leading at large wavelengths to an 
instability involving splay but not bend [T). One can 
likewise understand the results in Fig. Id for an exten- 
sile fluid, by sketching the expected nematic deforma- 
tion in front and at the rear of the particle (see inset in 
Fig. lb). Now the incipient instability is towards bend- 
ing, so pulling perpendicularly to the director gives the 
larger effect. Moreover, the active stress on the front 
of the particle now favors motion and the smaller one at 
the back opposes it: so drag is reduced. These arguments 
suggest that the v{F) curves should strongly depend on 
the anchoring conditions at the particle surface, and we 
have confirmed this numerically (data not shown). 

Building more quantitatively on these arguments, we 
now estimate the extra force acting on a moving particle 
arising from the active contribution to the stress (inte- 
grated over the colloidal surface). This depends on the 
surrounding flow field and nematic deformation, details 
of which can only be computed numerically. Nonetheless, 
we can argue that activity leads to an additional force 
on the particle which, on dimensional grounds, equals 
A{F)CR'^. To estimate A, we note that for F = the 
particle is immobile, so that the extra force should also 
vanish; A{0) = 0. To first order, dimensional analysis 
now suggests A ~ cF/K with c a positive dimensionless 
number of order unity [27]. (Our previous arguments 
show that A/F is positive for both signs of ("•) The force 
balance for the moving colloid in steady state is then 
given by F + c(FR^ /K = ^TrfjRv^ where r) is the passive 
nematic viscosity, which may depend on pulling direction. 
This equation together with Eq. [3] leads to the following 
prediction for the effective viscosity 
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FIG. 2: (Color online) (a) Plot of l/rjeE on i?^, yielding an 
approximately linear dependence in agreement with Eq. [2] 
(b) v{F) curve for a case with negative drag {R = 30, C = 
—0.002). Filled circles and open squares refer to simulations 
in which the force is applied to a quiescent and moving colloid 
respectively. The velocity fields in (c) and (d) correspond to 
the case of a quiescent particle (c) and of a particle moving 
to the left, opposite to the applied force (large arrow) (d). 



This provides a surprisingly robust explanation for the 
size dependence of the drag coefficient: Fig. 2a shows 
that for a colloid in a contractile fluid, plotting l/r^eff 
versus R^ yields a straight line as predicted by Eq. [4j By 
fitting our numerical data, we also obtain c ~ 0.06 for 
parallel pulling independent of which further validates 
our approximate theory. However Eq. [4] works less well 
in the extensile case. There the effect is smaller, and 
passive contributions to the stress tensor possibly play a 
more significant role. 

Our simplified theory in Eq. |4] formally predicts that 
when the dimensionless quantity cC^R? jK < — 1, 77eff < 0, 
giving a negative drag coefficient (^ = F/i' < 0). How- 
ever, one might reasonably expect that before this regime 
is reached our leading-order expression A ~ cF j K would 
break down and some different physics would come into 
play, so that ^ remains positive. Remarkably though, this 
is not the case. Although positive ^(F) appears to be re- 
stored at extremely small forces, for modest but finite 
F, we robustly find that the particle moves in the oppo- 
site direction to the externally applied force. A steady- 
state v(F) curve showing this effect for cC,R^ /K ~ 2.2 
is in Fig. 2b. This reverse-sigmoidal curve shows bista- 
bility at small forces, with two stable branches each of 
negative drag, ^ = F/v < 0, and also negative slope, 
dF/dV < 0. Interestingly, there seems to be no simple 
relation between negative local drag in contractile fluids 
and the homogeneous bulk rheology of active fluids. In 
the latter, a negative ratio of stress to strain rate can 
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arise, but only in the extensile case [3]. Moreover, the 
resulting flows are generically unstable [5]. 

To shed more light on the nature of the negative drag 
state in contractile fluids, we plot in Figs. 2c, d the active 
fluid flow around a stationary, unforced probe particle 
and that for one moving against the external force. For 
the static particle, a symmetric pattern of eddies is cre- 
ated by contractile forces arising from the deformation of 
the director (horizontal in far field), which splays around 
the sides of the probe and bends at its top and bottom. 
In the negative drag regime with the external force acting 
to the right, the eddies become asymmetric, creating a 
packet of left-moving fluid. The resulting leftward advec- 
tive velocity of the particle exceeds its rightward speed 
relative to the local fluid packet, giving an overall left- 
ward motion in the lab frame. 

We have examined how this remarkable steady-state 
flow pattern is reached dynamically. On introducing a 
rightward force, the particle first moves rightward (see 
Fig. 3) creating strong elastic distortions on its leading 
(right-hand) side. Such deformations then establish the 
leftward moving fluid packet (as in Fig. 2d), which sets 
up deformation in the surrounding director field to fi- 
nally create the steady state travelling flow packet with 
the particle, as well as the fluid, moving leftward. We em- 
phasize that this spontaneously moving state is quite dis- 
tinct from the bulk spontaneous flows known to occur in 
active fluids when activity exceeds a system-size depen- 
dent threshold [4 . Because we are below that threshold, 
in our simulations the fluid velocity is essentially local- 
ized in the neighborhood of the probe particle. 

A further exploration of the bistable negative-drag re- 
gion at small F is presented in Fig. 3. Here, we start 
with a static particle, of radius large enough to give neg- 
ative drag {c(R'^/K ^2.2) and first pull it to the right. 
The transient response was already described above, and 
leads finally to steady leftward motion. We now reduce 
F to zero: the particle keeps moving to the left, but slows 
down. We next start pulling the particle with a tiny left- 
ward force: the probe continues to move leftward, albeit 
even more slowly than before. If this last stage is re- 
peated with a much larger leftward force, the active flow 
and deformation pattern around the colloid dynamically 
reconstructs itself with the opposite sense, leading to a 
rightward probe velocity and once more to a negative 
drag. This sequence of states can be viewed as a tra- 
jectory on the v{F) curve of Fig. 2b. The initial large 
rightward force puts the particle on the lower branch of 
the curve to the right of the vertical axis. As the force is 
removed and then applied to the left, the system ascends 
this lower branch. Where that branch ends, the system 
must jump to the upper branch. 

Conclusion: we have simulated numerically a simple 
microrheological experiment, which should be realizable 
in the laboratory (with some caveats, noted below). In 
this experiment a colloidal probe is pulled through an 
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FIG. 3: (Color online.) Dependence of the external force 
(step functions, axis on the right) and of the particle velocity 
(axis on the left) on time for two different 'experiments'. In 
both cases a force F = 0.8 is initially applied to the particle 
and turned off. Then a smaller force is applied along the 
direction of particle motion: the blue dotted line and the red 
dot-dashed one give the values of the external force in the two 
cases. The correspondent particle velocities are given by the 
black solid and green dashed lines respectively. The dynamics 
is discussed in the text. In the sketches arrows on the fluid 
refer to the active force direction, the arrow on the colloid 
represents the external force, while that above the particle 
shows the actual direction of motion. 

active nematic fluid. We have shown that the activity 
leads to a non-Stokesian drag force that increases approx- 
imately quadratically with particle size. This behavior 
creates a new regime, arising in a contractile fluid such as 
an actomyosin gel, at large values of a dimensionless pa- 
rameter c(R^ /K. Strikingly, in this regime the colloidal 
probe is predicted to move against the driving force to 
create a stable steady state of negative drag. This con- 
trasts with the bulk behavior of active fluids where neg- 
ative ratios of stress to strain rate can arise in principle, 
but are unstable, and expected only for the opposite sign 
of activity (extensile rather than contractile) [3 . 

Stable negative drag is counterintuitive but no physi- 
cal law prevents it in active systems. Related phenomena 
have been reported for an ensemble of molecular motors 
with load-accelerated dissociation ^8], in axon mechan- 
ics [29 , for filament fluctuations in active media [7 , and 
also in the upstream migration against a flow field of 
slime- mold cells [W and bacteria [31]. In our context, 
the anomalous (ultimately negative) drag stems from the 
active stress, Eq. |2j that emerges from a well-accepted 
coarse-grained description of activity in systems such as 
cytoskeletal gels [1]. Finally, an important requirement 
for any experimental test of our predictions with acto- 
myosin is that issues arising from macroscopic network 
clustering should be avoided [32]. This might be feasible 
either through a suitable choice of parameters such as 
myosin concentration etc., or by working in metastable 
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uniform networks. However, such experiments may well 
prove to be challenging. 
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